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Abstract 

We discuss the phenomenology of neutrino mixing with bulk fermions in the 
q \ context of supernova physics. The constraints on the parameter space fol- 

CNl ■ lowing from the usual energy loss argument can be relaxed by four orders of 

magnitude due to a feedback mechanism that takes place in a broad region of 
the parameter space. Such a mechanism also affects the protoneutron star evo- 
lution through a non trivial interplay with neutrino diffusion. The consistency 
with the SN 1987A signal is discussed, as well as the implications for delep- 
tonization, cooling, composition of the neutrino flux and the delayed explosion 
scenario. 
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^ ■ 1 Introduction 

O 

Sterile neutrinos from extra dimensions are not likely to play a primary role in atmospheric 
and solar neutrinos fl:. On the contrary, astrophysics, cosmology and rare decays still 
represent natural stages for their exotic performances. In this paper, we focus on the 
implications of a possible mixing of the electron neutrino with a Kaluza-Klein (KK) tower 
i-Eh ■ of sterile neutrinos in supernova (SN) physics. 

The potential relevance of bulk fermions for neutrino physics has first been pointed 
out and analyzed in [21 El- Fermions propagating in large extra dimensions appear as 
light, dense KK towers of sterile neutrinos in 4 dimensions, possibly mixing with the 
Standard Model ones. In turn, the existence of extra dimensions much larger than the 
inverse electroweak scale is natural in string inspired models with a fundamental scale in 
the TeV range pEJE], if the standard relation between that scale and the Planck scale is 
implemented. 

The basic ideas about supernova explosion and cooling, nicely confirmed by the neu- 
trino signal from SN 1987A [B], translate into significant constraints on the energy loss 
through invisible channels [Jj. In fact, in order for the observed neutrino signal to be 
accounted for, the invisible energy loss should not cool the protoneutron star on a time 
scale shorter than about 10 seconds. The bulk of a large extra dimension provides at 
least two interesting examples of invisible channel: KK graviton emission j^j and mixing 
of Standard Model and bulk neutrinos, the possibility we are interested in. With respect 
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to the case of a single sterile neutrino [S], conversion into KK towers of sterile neutrinos 
is enhanced by the higher number of available states and especially by matter effects [10] . 
The constraints that follow on the parameter space of bulk neutrino models look at first 
sight quite severe ^U]. However, those constraints are relaxed by an interesting feedback 
mechanism 1 that prevents unacceptable energy loss ^2] and affects the protoneutron star 
deleptonization and cooling through a non trivial interplay with diffusion. By suppress- 
ing the neutrino MSW potential and then making it negative everywhere, the mechanism 
leads to a self-limiting of the potentially dangerous, MSW-enhanced, energy loss into the 
bulk and requires a revisitation of the bounds in the literature. The consequences for SN 
physics, in particular explosion, deleptonization, cooling and spectra, are not less inter- 
esting. From this point of view what we discuss here is an example of how the standard 
ideas about the protoneutron star evolution can be affected in presence of new physics in 
an unexpected and intriguing way. 

The paper is structured as follows. In Section [2] we shortly review the salient points 
about (electron) neutrino oscillations into bulk neutrinos. Section |3] is the central part of 
the paper. Since a reliable study requires taking into account diffusion, we first set up in 
Section fc.ll a toy model of the protoneutron star core which, although minimal, contains all 
the neutrino transport physics relevant to our discussion. In Section \'A. 21 we incorporate in 
the model the effect of electron neutrino and antineutrino conversion into bulk neutrinos. 
The effect on the evolution of the protoneutron star core is described in Section 13.31 The 
consistency with SN 1987A is shown in Section 13.41 where the implications for the neu- 
trino luminosity, the flux spectrum and the shock reheating are also discussed. Section 
summarizes the constraints on sub leading non self-limiting conversion probabilities, that 
provide one of the relevant bounds on our parameter space. We conclude in Section [3 

2 Framework 
2.1 The supernova 

Before embarking in the detailed analysis, let us shortly and qualitatively review the 
standard ideas about SN explosion and cooling [TUl HI] . 

Type II Supernovae originate from evolved, massive stars whose iron core reaches 
the Chandrasekhar limit and collapses under gravitational pressure. During the collapse, 
neutrinos start to be produced by electron capture and get trapped when densities of 
order 10 12 gem -3 are reached. The inner part of the iron core has an homologous collapse 
that stops when, after a fraction of a second, stiff nuclear densities of order po = 3 • 
10 14 gcm~ 3 are reached. The free falling outer material then builds up a shock that 
propagates outwards. It is believed that this "prompt" shock has not enough energy to 
give rise to a successful explosion and hence stalls after ~ 0.1 sec at a radius of order 
100 km. Below the shock, the protoneutron star has a high density inner part, the inner 
core, with a radius of 10-30 km and a bloated outer part, the mantle, that accretes mass 
falling through the shock. Let us leave for a moment the shock to its uncertain fate and 
follow the evolution of the protoneutron star below it. Most of the gravitational binding 
energy of the progenitor star, a few 10 53 erg, is stored there. Only O (1%) out of that will 

1 A different type of feedback was considered in 



2 



go into explosion kinetic energy or electromagnetic radiation. Most of it will be released 
through neutrino emission in a few tens of seconds. 

The stiff inner core has a mass of order of the solar mass, M Q ~ 2 • 10 33 g, a density 
of order po, a temperature of 10-30 MeV and a lepton fraction per baryon Yl ~ 0.35. 
Electrons and electron neutrinos are highly degenerate in such conditions, with Fermi 
momenta of order 300 and 200 MeV respectively. The mean free path of a 200 MeV 
neutrino is A ~ 10 cm. This determines the time scale for neutrino diffusion and emission: 
^diff ~ 3i? 2 /A ~ 10 sec for a neutrino at a depth R = 10 km in the inner core. Analogously, 
the total neutrino flux and spectrum can be roughly estimated from the knowledge of the 
star structure and temperature where neutrinos experience their last energy exchange. The 
mantle looses most of its lepton number, accretes most of its mass, cools and contracts in 
only 0.5-1 sec. Then, as it cools, the whole protoneutron star slowly contracts to a radius 
of order 10 km. After a few tens of seconds, we are left with a proper neutron star. 

Since supernovae do explode, we can be confident that the shock has in the meanwhile 
successfully expelled the outer matter. An important contribution to the revival of the 
stalled shock is given by the energy deposition of the neutrino flow. A O (1%) fraction 
of the energy flowing might be enough to give rise to a sufficiently energetic explosion 
a second or so after collapse. Critical are of course the efficiency of the energy transfer 
mechanism and the intensity of the neutrino source. 

Although this "delayed shock" mechanism is currently believed to play an important 
role, reproducing the explosion in computer simulations it is still non trivial at present |15j . 
This does not necessarily mean that something is missing in our general understanding. 
The problem is likely to originate from the complexity of the system to be modeled and 
the variety of the physics involved. It may just be necessary to turn to full 3D simulations 
able to account for non-spherically symmetric effects. Convection is likely to be important 
in this context, since it may boost the neutrino luminosity 2 . Moreover, the neutrino 
energy deposition in what was the prompt shock is likely to be just one of the several 
effects giving rise to the explosion |17j and the outcome critically depends on a number 
of poorly known variables like the progenitor structure, the equation of state, the details 
of neutrino transport, etc. We also note that a purely neutrino driven explosion is not 
obviously reconciliable with the indications of a non-spherical mechanism |18| 3 . Still, it 
is worth pointing out that the new physics and mechanism we discuss in this paper might 
also play an important role by enhancing the neutrino luminosity. This will be discussed 
in Section 

The basic ideas about protoneutron star deleptonization and cooling have been nicely 
confirmed by the observation of a neutrino signal in the Superkamiokande and 1MB de- 
tectors in coincidence with SN 1987A [§]. Despite the low statistics, there is definitely a 
significant agreement between the measured neutrino flux duration and intensity and the 
estimates above. This allows to set limits on drastic departures from the described picture, 
in particular on invisible energy loss channels that would waste the available energy on 
timescales shorter than 10 sec. 

2 However, whether such effects do provide a solution of the SN problem is still an open issue |1^J. 
3 The invisible channel we are discussing might itself provide a source of asymmetry. 
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2.2 Bulk neutrinos and feedback 



We consider a Kaluza-Klein (KK) tower of sterile neutrinos living in the bulk of a single, 
large, flat, gravitational extra dimension. Additional smaller dimensions will be in general 
present. What we consider is therefore an effective theory valid up to scales below the 
inverse scales of the smaller dimensions. The physics we are interested in involves bulk 
neutrinos with masses within ~ 100 keV. We therefore assume that the inverse scale of 
the additional smaller dimension is larger than 100 keV and neglect them in the following. 

As for the size R of the extra dimension we consider, we only require it to be in the 
range 

10- 3 eV<i<lkeV. (1) 

The lower bound comes from tests of the Newton's law at small distances 4 and the origin 
of the upper bound, as well as the bounds below, will become clear by the end of the 
section. 

Apart from R, which gives their density of states, all we need to know about bulk 
neutrinos is their mixing with the electron neutrino. We assume that the mixing Ok of v e 
with the (almost) sterile fc-th mass eigenstate N\. is inversely proportional to its mass M^, 

e k ~ , ( 2 ) 

a generic feature of most models in the literature EH C2 EH EOJ • Eq. © defines the 
mixing parameter m. 5 

We are interested in a range for m such that mR <C 1, which ensures that all mixings 
6k are small, and 

10~ 12 eV < m 2 R < 10" 4 eV . (3) 

The parameter m 2 R sets the scale of all the effects we discuss. The upper bound ensures 
the smallness of the transition probabilities over a mean free path and values of m 2 R 
smaller than the lower bound give unobservable effects. The range in eq. © is further 
restricted by the expectation that subleading contributions to the oscillation probability 
give rise to the bound mR < 10 -5 , as we will discuss in Section |IJ Together with the 
upper bound on 1/R, the latter gives m 2 R < 10 -8 eV. We will therefore concentrate on 
the first four orders of magnitude of the range © • 

For appropriate sizes of the volume associated to the smaller dimensions, all values of 
1/R, m 2 R in the ranges above are consistent with the relation (Mpi/M*) 2 = FJ^ 2-nRiM* 
between the volume of the extra space with the 4D Planck mass, for appropriate sizes of the 
smaller dimensions. At the same time, those values are compatible with cosmological and 

4 Depending on the specific model and on the scale of the mixing term m in eq. (|5J below, it may be 
necessary to strengthen the lower bound to 1/72 > 10 -1 eV in order to avoid undesired effects in solar 
and atmospheric neutrino experiments. In fact, for 1/R < 10 _1 eV, the first modes of the KK tower may 

participate to oscillations if their masses are given by n/R, n = 0, 1, 2, This is not the case in presence 

of a Dirac or Majorana bulk mass fi, since in this case all bulk neutrinos are heavier than \x [191112) . Values 
of 1/R as low as 10~ 3 eV would then be automatically compatible with not having effects in the observed 
oscillation signals. Provided that fi <C keV, such a mass term would not affect SN physics. 

5 The exact dependence of Mk on 1/R may depend on the specific model. We only assume the density 
of states dk/dM to be R at scales around lOkeV. The parameter m usually originates from Yukawa 
interactions involving the 5D neutrino, the lepton doublet and an Higgs field. 
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astrophysical bounds from KK graviton production before Big-Bang nucleosyntesis |21| 
and in supernovae |Hj. 

Having described the neutrino parameter space, we now discuss neutrino oscillations 
in the inner core of the protoneutron star. Oscillations are dominated by matter effects. 
While bulk neutrinos are sterile and therefore unaffected by matter, the electron neutrinos 
acquire an effective squared mass = 2EV, where 



V = V2G F n B Qy e + 2Y Ue - M 



(4) 



is the MSW matter potential, ns is the baryon number density and Y e<Ue = n ejVe /riB is 
the number fraction of e, u e per baryon. More precisely, n ejVe = N e - Ve — N e + Pe denotes 
the difference of the particle and antiparticle number densities. Also, we denote by Yl the 
total lepton fraction Y e + Y Ve . 

The neutrino contribution to V is important in the inner core. Also important is the 
negative contribution from neutrons, that allows the potential to have both signs, as shown 
in Fig. ^ According to the sign of the potential, matter effects will enhance neutrino or 
antineutrino disappearance in the bulk. 

The MSW potential ranges between — 10 eV and 5eV. Correspondingly, in the early 
stages of evolution, most electron neutrinos have an effective mass in the range 

(lOkeV) 2 — I— < m 2 cS < (50keV) 2 L 



100 MeV ~ ctt ~ v '100 MeV ' 

where E is the neutrino energy, as shown in fig. A resonance is met when the effective 
mass matches the mass of a state in the KK tower. The condition mR <C 1 ensures that 
the width of the resonance is small compared to the distance between levels. 

When an electron neutrino emerges from an interaction in the inner core, its mixing 
with the massive states in the KK tower depends on its effective mass m e g and therefore 
on the values of ng, Y e , Y Ue where the interaction took place. If V is positive, v e is 
predominantly made of the mass eigenstate Nj. such that < m c s < M^+i, with only 
a small component from the other eigenstates 6 . If V < 0, on the other hand, v e does not 
mix with the sterile states — the mixing takes place in the antineutrino sector in this case. 
After its interaction, the neutrino travels a distance of order of its mean free path, say 
about 10 cm, before experiencing the next interaction. Let us determine the probability 
P{y e — ► bulk) that the electron neutrino has in the meanwhile oscillated into a bulk 
neutrino assuming that V is positive. If V is negative, the probability is negligible and what 
below applies to antineutrinos. Over 10cm, the values of ub, Y e , Y Ve change only slightly, 
but the corresponding change in the effective mass may be sufficient for the neutrino to 
predominantly mix with a different mass eigenstate Nh- The dominant contribution to 
the disappearance probability P{y e — > bulk) is then given by the probability that Nf. has 
not turned into Nh. This level "non-crossing" probability can be computed by means of 
the Landau-Zener formula [5] 

P{y e -► bulk) ~ 1 - e"5 |An|7 ~ -| An| 7 , (5) 



6 This is always true except when m e g coincides with some Mk within a resonance width. In such a 
case, u e will predominantly be made by a superposition of Nk and iV^+x- See |12| for a detailed discussion 
of mixing. 
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where An is the number of levels crossed, 7 is the adiabaticity parameter and we assumed 
P(y e — > bulk) <C 1, as follows from m 2 R <C 10 eV (for P(u e — > bulk) ~ 1, the analysis 
is complicated by non linear effects and back conversion of bulk neutrinos into u e ). The 
Landau-Zener formula can be extended to the case of multiple resonances because the 
latter are well separated. 

The adiabaticity parameter is given by 

7 = ^ ( 6) 

1 E{dV/dr) cos <j> ' K ' 

where <j) is the angle between the neutrino momentum and the radial direction. The 
number of levels crossed over a distance L is 7 

dV RE T 

An = — Lcoscp. (7) 



dr m e ff 



We therefore conclude that 



1/2 



P{v e -+ bulk) ~ 6{V)L ^rn'R 1^1 (8a) 

/ lyl \ 1/2 

P(^^bulk)^0(-y)L^=m 2 i?(^yj . (8b) 

The step function accounts for the fact that neutrinos (antineutrinos) do not cross any 
resonance if V < (V > 0). Correspondingly, we have neutrino disappearance for V > 
and antineutrino disappearance for V < 0. The probabilities depend on the neutrino 
parameters only through m 2 R, to which they are proportional. The upper bound in 
eq. © ensures that P(u e — ► bulk) -C 1. 

In obtaining the previous equation, we have approximated the electron neutrino with 
the mass eigenstate of which it is predominantly made. Additional subleading contribu- 
tions to P{y e — > bulk) come from the small mixing with the neglected eigenstates, that 
amount to a fraction of order mR of the electron neutrino. Their effect will be discussed 
in Section 2] 

Let us now discuss the motivation for restricting the range of 1/R to l/R < keV. 
The number of resonances crossed in a mean free path can be estimated from eq. ( |7| ). 
Roughly speaking, for 1/R < eV most neutrinos are likely to cross at least a resonance 
in a mean free path — hundreds in the lower part of the allowed range. For 1/R> eV, 
on the other hand, resonances are met only by neutrinos with energies within ±5E from 
a set of discrete values E^ = M%/(2V) in which m e g = The discrete values E^ are 
separated by AE k = M k /(RV). When 1/R > eV, SE = ELcos<p{dV/dr)/V is smaller 
than the separation AE k and the resonances are not crossed for arbitrary values of the 
energy anymore. However, there will always be neutrinos crossing a resonance as long as 
some of the energies E k fall in the neutrino spectrum. The condition 1/R < 1 keV ensures 
that there are at least 20-50 points in the spectrum where neutrinos cross a resonance 
at the beginning of the evolution. It also ensures that there will be a few wells in the 
spectrum even when, at later stages, the product V ■ E will decrease by a factor 1000. 

7 For 1/7? > 1 eV, not all neutrinos cross a resonance and An should be considered as the average 
number of level crossed (see below) . 
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We now have all the ingredients to understand how v e conversion in the inner core 
tends to draw the MSW potential to zero everywhere through a feedback mechanism. Let 
us consider a region in the inner core where the MSW potential is positive. Eq. (JHJ) shows 
that the neutrino conversion probability is proportional to VV. On the other hand, once 
a neutrino has been lost in the bulk, V decreases. This is because V oc 3Y e /2 + 2Y Ue — 1/2. 
As a consequence of the loss, the electron and electron neutrino number densities will 
readjust to restore (5 equilibrium, but 3Y e /2 + 2Y Ue will certainly decrease. The more 
neutrinos are lost in the bulk the smaller V becomes. It may happen, and we will see that 
it does happen, that V approaches zero before all the lepton number and energy are lost. 
At this point the neutrino conversion through the MSW enhanced probability (jS^) stops. 
In the region where the potential is negative, antineutrino escape would take place (at a 
much slower rate) drawing again the potential toward zero. As we will see, the possibility 
of reaching zero depends in this case on the local values of the thermodynamic variables. 

The effect we are dealing with is potentially quite interesting. Besides limiting the 
energy loss, in fact, the condition V = represents an important constraint on the relative 
abundances, with significant implications for the SN phenomenology. However, several 
issues have to be addressed. First of all, the V = condition is spoiled by neutrino 
diffusion, which tends to lower V below zero. Antineutrino conversion into bulk neutrinos 
then tends to bring V back to zero at the expenses of an additional energy loss. The actual 
evolution will be determined by an interplay of diffusion and conversion. Moreover, it has 
to be proven that the mechanism is compatible with SN 1987A observation. Does the inner 
core loose all its lepton number and energy in the attempt of reaching the V = condition? 
Or while neutrinos try to diffuse out? How quickly is the condition approached? Is the 
escape of the rare antineutrinos efficient enough? Qualitative arguments are not sufficient 
to reliably address the issues above, not to mention attacking the problem of implications 
for possible future SN observations. On the other hand, a full, detailed simulation of the 
protoneutron star evolution is clearly beyond the scope of this paper. We therefore set up 
in the next section a toy model of inner core evolution which only incorporates the basic 
ingredients of neutrino diffusion but contains all the physics relevant to our study. 

3 Inner core evolution 

3.1 Standard evolution: equations and assumptions 

We now discuss the equations we use to model the evolution of the inner core. We start 
from the case of no mixing with bulk neutrinos. We concentrate on the inner core because 
the bulk of energy and lepton number is stored there. Moreover, the loss rates and the 
effectiveness of the feedback on the MSW potential are highest in the core. That is because 
the neutrino densities are large and because the transition probability in eqs. (jHJ grows 
with the potential, which is proportional to the baryon density 8 . 

Focusing on the inner core has also practical advantages. Unlike the mantle, the 
core settles into local thermodynamic equilibrium very quickly after the collapse. The 
hydrodynamics is also much simpler, with the mass profile becoming essentially constant 

8 The transition probability goes with (V/E) 1 ' 2 , where E also grows with the baryon density ub- 

1/3 

However, V <x tib, whereas E oc n fl . 
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in a few hundreds milliseconds ^3]. The mantle of the protoneutron, on the other hand, 
accretes matter for the first 0.5-1 sec and then slowly contracts. 

Assuming thermodynamic equilibrium and neglecting general relativity effects, the 
basic equations of neutrino transport in absence of mixing are [221 

V • F L (9a) 

V • F e + /V? • F L , (9b) 

where s is the entropy per baryon and Fl = F u — Fp, F 6L = F 6u + F e - are the lepton 
number and lepton energy density currents respectively. The neutrino currents are 

Fu = ~\ £ f-^K(E)VUE) (10a) 

p e v = -\ E / E Kw&um , (iob) 

where p is the neutrino momentum, E = p is the neutrino energy, X u is the neutrino 
mean free path (mfp), which also depends on the local thermodynamic variables, and 
f v is the Fermi-Dirac distribution, f v {E) = ( y e^ E ~^/ T + l) -1 , which depends on the 
temperature T and the neutrino chemical potential \i v . Analogous expressions hold for 
antineutrinos with \i v — * —\l v . We assume that muon and tau neutrinos have vanishing 
chemical potential 9 . As a consequence, they do not contribute to the lepton number 
current and only give a thermal contribution to the energy current. The neutrino chemical 
potential, as the electron, neutron and proton ones, p e , \i n and be obtained in 

terms of the thermodynamic variables T, Yl and p (p is the mass density), by solving the 
equilibrium equation p e — p u = p n — p p . Since we do not solve the whole protoneutron star 
evolution, we also have to specify boundary conditions for eqs. ©. We approximate them 
by imposing that the lepton number and energy fluxes at the border are proportional to 
the neutrino number and energy density. We assume spherical symmetry. 

The transport equations © for Yl and T incorporate all the physics we are interested 
in before inclusion of non-standard effects. Their solution however requires i) the knowl- 
edge of the mean free paths, ii) the knowledge of the equation of state, which enters the 
equilibrium condition and the entropy, and iii) the determination of the density p(r, t). In 
order to prune our considerations from unnecessary complications and concentrate on the 
transport physics of eqs. @, we address the three points above in the simplest possible 
way. First of all, we take advantage of the hydrodynamic stability of the inner core by 
using a constant (in time) density profile p(r). As for the equation of state, we consider 
a core made only of n, p, e^, v e ,ti,T, ^,/i,t, 7 and we take into account the strong nuclear 
interactions through a modified dispersion relation for the nucleons. That is, the nucleons 
are approximated as a perfect Fermi gas, with an effective mass m* N carrying a depen- 
dence on the state variables. A detailed calculation of m* N has to rely on some effective 
theory of nucleon-nucleon interaction and, in general, will deviate substantially from the 

9 Given the high energies involved, a non- vanishing number of muons could participate to the inner core 
life, giving a non-vanishing chemical potential. 
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vacuum value = 939 MeV. The result of such theories can be reproduced by a simple 
expression carrying only a dependence on the medium density JH] 

„ _ m N 
mN -l + f3 0P / P[ 

where @q is chosen to be 0.5 and po = 3 10 14 gr/cm 3 is the reference nuclear density. 

Finally, we have to define a prescription for the mean free path of neutrinos of different 
flavors. For electron neutrinos the most important contributions to the opacity come from 
the absorption reaction v e +n — * e~ +p and from the scattering on neutrons. Scattering on 
protons and, to a smaller extent, on electrons, are also sizable. For za/s and f T 's, only the 
scattering reactions are important. A detailed analysis of neutrino opacities should take 
into proper account the cross sections for these and several other subleading processes, as 
well as the status of the nuclear matter (degenerate or non-degenerate nucleons, consequent 
Fermi blocking effects), the nucleon-nucleon interactions, the multicomponent nature of 
dense matter (e.g. the presence of hyperons), etc. For a comprehensive discussion on this 
subject we refer the reader to [SB] an d references therein. For our purposes, it is sufficient 
to model the neutrino mean free path with a simple inverse quadratic dependence on the 
neutrino energy 

E 2 E 2 
X Va (E) = X° Ue -^ , A„ M) „ r (£) = A^ ^-^l , (12) 

where we take A° = 1.2 cm and A" v = 4.8 cm at the reference energy Eq = 260 MeV. 
The values of A^ e and A° Vt determine the diffusion time scale. When compared with the 
exact general expression of ref. J23, (taken at typical density, temperature and leptonic 
fraction), the above proves to be a good approximation for the upper portion of the range 
of neutrinos' energies of our interest while it turns out to be poorer in the lower part. 
Moreover, the above choices yield an evolution whose timescale agrees with the results of 
more sophisticated analysis |24j . 

With the choice above on the energy dependence of the mean free paths and setting 
= A y in the small antineutrino contribution, one finds the simple expressions 

2 

-F L = a e Vp u , -4 = yV^ + (a e + a M + a T )yVT 2 , (13) 

where Oj = \® Eq / (6ir 2 ) , i = e,fj,,r. 

As mentioned, the simplified model we use for discussing the inner core evolution only 
contains the basic features of neutrino diffusion relevant to our discussion. It is certainly 
not suitable for addressing more complex issues such as non-spherically symmetric effects 
likely to be involved in the explosion. 

3.2 Neutrino escape into the bulk 

Eqs. ([9113ft . together with the equilibrium equations and the boundary conditions represent 
our simplified model of inner core deleptonization and cooling in absence of non standard 
dynamics. We now incorporate the effect of neutrino and antineutrino escape in the bulk. 

Electron neutrinos with energy E contribute to the lepton number loss rate by (P) /A 
each, where P = P(y e — ► bulk) and "()" indicates an average over the distance traveled 

9 



(11) 



L. 10 Since in our case the probability is linear in L, we simply have (P) /A = P/L = 
6{V)tt / {2^)m 2 R(\V\/ E) 1 / 2 . The total neutrino number loss rate 1^ follows fr om inte- 
gration over neutrino momenta and can be written by means of the Fermi integrals F a 

as 

>(V)_2„ /iT77^5/2„ „ 9(10 



OC 



Analogously for antineutrinos. Putting all together, we find the following modified evolu- 
tion equations for Yl and T: 

HB ^r = ^ " {ae ^ ~ i^ m2R ^ 1 t5/2jFs/2 (15a) 



(15b) 



where cr = sign(V). We recognize in eq. (|15b) the heating term a e (V^,) 2 associated to the 
degradation of the degeneracy energy of the neutrinos reaching regions with lower chemical 
potential. This term partially originates from the chemical potential contribution to T ds. 
With the same origin, there is a heating term in the energy loss part proportional to F%j2 
that can overcome the cooling term proportional to F 5 / 2 - I n this case the temperature 
increases because of the degradation of the degeneracy energy of neutrinos at the Fermi 
surface that downscatter to replace neutrinos inside the Fermi sphere escaped in the bulk. 

In eqs. (fT5]l the potential V, as p v and s, is a function of the thermodynamic variables 
p, T, Yl. A simple approximation is useful to understand the dependence of V on those 
variables, in particular Yj,. By neglecting the fi n — pp term in the beta equilibrium equation 
Pe — fJ-u = Pn — Mp> one obtains Y e = 2/3 Yl, Y u = 1/3 Yl- Such a rough approximation 
turns out to be pretty good for V, where errors cancel to give, with an accuracy of 5% or 
less, 

V ~ V2G F n B 0y L - , (16) 

which vanishes for Y® ~ 3/10 independently of T, p. Eq. ()16|) illustrates what happens 
to the MSW potential when a neutrino or an antineutrino is lost in the unit volume. The 
lepton number lost or gained is redistributed by the beta equilibrium between neutrinos 
and electrons and the net effect is a variation of the potential by the amount :f (5/3)y2Gp. 
Moreover, through eq. (|16|) it is possible to explicitly see the effect of the feedback. Let us 
consider the limit in which the loss term dominates the evolution and diffusion is negligible. 
Then, when Y L > 3/10, the RHS is negative and Y L decreases. When Y L < 3/10 the RHS 
is positive and Yl increases. When Yl = 3/10, the RHS vanishes and the lepton fraction is 

10 Note the difference between (P) / A and {P/L}. While the latter expression looks at first sight more 
correct, it actually does not take into account, as the former does, the variation of the number densities 
in the time necessary to travel the distance L. 
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constant. In the approximation above, 3/10 is the value of Yl at which the non standard 
terms in eqs. (|15[1 aim. We will see that for m 2 R 3> 10 -12 eV this value is indeed reached 
and maintained for some time in some region of the inner core. In other regions, the value 
is not reached but its attraction on the lepton fraction affects the deleptonization. In our 
numerical implementation of eqs. (|15jl. we do not make use of the approximation (|lri[). The 

V = condition then becomes a constraint on p, T, Yl that translates into Yl = Yj*(p, T). 

3.3 Phases of evolution 

We start to follow the evolution of the inner core a few hundred milliseconds after collapse. 
We consider a core of 1.5 M & with a radius of 12.7km and a typical density profile described 
by p(r) = p c /(l + (r/f) 3 ), with p c = 7.5 10 14 gr/cm 3 , 4/3 7r f 3 p c = 1.1 M . The initial 
profiles for T and Yl are shown in fig. (the thick dashed lines). Their main features follow 
from models of core collapse |24j . The detailed structure of those profiles is not important 
for our purposes. Diffusion will in fact soon smooth them. Moreover, the mixing with 
bulk neutrinos, when important, also affects the early stages of SN evolution. This point 
will be discussed in more detail below. 

Three processes contribute to the variation of the lepton fraction: diffusion, neutrino 
and antineutrino conversion into bulk neutrinos. Each of them has a characteristic time 
scale, tdifib tu and t v respectively, defined e.g. as the elapsed time per fraction of Yl 
variation. While t u turns out to be finite and roughly the same in all the region where 

V > 0, tp strongly depends on the position in the V < region and may become infinite, 
as we will see. In any case, t v <C tp. In most of the inner core, in fact, the beta equilibrium 
forces p u to be positive and the neutrinos to be degenerate, p u 3> T. As a consequence, 
the antineutrino number is strongly suppressed and so is the lepton number variation due 
to antineutrino escape. 

While the time scale for diffusion is typically ~ 10 sec, the escape time scale depends 
on m 2 R. Moreover, neglecting diffusion, the evolution due to the escape depends on m?R 
only through an overall scale factor, Y^(r, t; k m 2 R) = Yl(v, kt;m 2 R), T(r,t; k m 2 R) = 
T(r, kt; m 2 R). The larger m 2 R, the faster the variation of Y^: 11 



We therefore distinguish two regimes: 

• m 2 R <C 10~ 12 eV, or t„ 3> t^g. Neutrino and antineutrino escape are too slow to 
affect significantly the evolution of the protoneutron star. Consistency with the 
observed neutrino signal from SN 1987A is ensured. 

• m 2 R 3> 10~ 12 eV, or t v <C tdifr- Neutrino escape changes Yl faster than diffusion. 
The consistency with the SN 1987A signal would be in question if the feedback were 
not taken into account. Deleptonization and cooling are affected. 

In the remaining of this subsection, we will concentrate on the case of fast escape, say 



n For m 2 R > 10 7 the thermal processes would not be fast enough to maintain equilibrium. The time 
scale above would be affected, but the results below would not. 




(17) 



2 R > 10" 11 eV. 
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We can distinguish two phases in the evolution. In the first one neutrino escape dom- 
inates and diffusion can be neglected. The second is a mixed phase where antineutrino 
escape dominates in some region and for some time and is comparable to diffusion else- 
where. 

The first phase takes place in the V > region, corresponding to Yl > Y?, on the 
left of the vertical lines in fig. |21 The change of Yl and T due to antineutrino conversion 
(where V < 0) and diffusion is negligible, since the time scales are slower. Due to neutrino 
disappearance, both Yl and the energy density start decreasing. The essential point is 
whether the value of Yl that stops the MSW conversion is reached before all the energy 
has been lost. This turns out to be the case in all the region where V > 0. In fact, by 
neglecting the diffusion terms in eqs. (|15jl we find the following relation between lepton 
fraction and the entropy variation: 

ds = ( F b / 2 {^ u /T) /^A _2/^, 

0Y L ~ \F 3/2 (^/T) TJ ~ 7T ' 1 j 

where the partial derivative is taken at a fixed position and the approximation Os/OYl ~ 
— 2/7 Hv/T holds for degenerate neutrinos. Since such an approximation does hold for 
Yl > 0.3 and T > 30MeV, eq. ((TH|) shows that the temperature in the V > region 
actually grows while it deleptonizes. In turn, that means that the V = condition is 
comfortably attained spending only a fraction of the available energy. Within a time 
t = 0.5-1 sec (10~ n eV/m 2 R), Y L has reached Y£ in all the V > region and the fastest 
phase of evolution has ended. The profiles are now shown as solid thick lines in fig. |HI The 
temperature has increased as a consequence of the conversion of degeneracy energy into 
thermal energy. As mentioned in Section 13.21 such a conversion takes place when states 
deep inside the neutrino Fermi sphere are emptied by neutrinos escaping in the bulk. In 
summary, after the initially dramatic energy leak into the bulk, most of the energy in the 
V > region gets locked as Yl reaches Y®, which happens on the t v time scale. The 
thick lines in fig. 01 can therefore be considered as the initial condition for the subsequent 
evolution. 

Before discussing the next phase, a comment on the initial condition we used is in 
order. As mentioned at the beginning of the subsection, for m 2 R 10 _12 eV the initial 
profiles are likely to be affected by the neutrino escape. For m?R = 10 -8 eV, for example, 
moving from the initial to the Yl = Y® profile only takes about a millisecond. Therefore, 
the initial profiles in fig. |HJ following from analysis that do not take into account the new 
effect, cannot be trusted hundreds milliseconds after core bounce. On the other hand, 
given the fixed point character of the evolution in the fast phase, we can be confident 
that the thick solid profile is indeed reached as soon as the core settles. We are therefore 
entitled to use that profile in the subsequent evolution. Why not to use it in the first place 
and ignore the fast phase of evolution, then? Because the simulation based on eq. (|18|). 
although it does not take into account the conditions met in the early stages of the SN, still 
gives a conservative estimate of the temperature variation and of the amount of the energy 
loss, a key piece of information for what follows. In summary, the thick profiles in fig. |3] 
can be considered realistic initial profiles for the next phase within the approximations 
and uncertainties associated to our approach and to core collapse models. Needless to 
say, the original initial profiles are suitable as initial conditions for the standard diffusive 
regime (m 2 R < 10~ 12 case) within the same uncertainties. 
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Figure 1: Profile of the MSW potential 
of electron neutrinos corresponding to the 
initial profiles in fig. 



Figure 2: Effective mass of a 100 MeV 
neutrino (where V > 0) or antineutrino 
(where V < 0) induced by the potential in 

fig.m 



After neutrinos in what was initially the V > region have been locked, antineutrinos 
loss begins to be significant in the V < region. This happens on the slower time scale 
tp, which depends on the position in the core. Since antineutrinos are lost, the lepton 
fraction Yl, initially smaller than Y®, grows and V becomes less and less negative. Again, 
the essential point is whether the V = condition, or Yl = Y®, is reached before all the 
available energy is lost. The relevant equation is in this case 

ds _ (F 5/2 (-^/T) + M /a* + 6\ (19) 



9Y L \F V2 {-ix u /T) TJ \T 2 

The loss turns out to be more pronounced than in the previous phase. Actually, the energy 
loss per unit of lepton number change is smaller, since the average antineutrino energy is 
5/2 T whereas the average neutrino energy was 5/7fi u . However, the entropy variation is 
larger since the increase of Yl leads to an increase of the energy stored in the degenerate 
neutrino or electron sea at the expenses of the temperature and entropy. That is why 
the RHS in eq. (j!9[) is larger than in eq. (|18|). Furthermore, the change in Yl needed to 
reach Y® is larger in the external part of the inner core, where Yl is lower. In order to 
clarify the situation, we need to numerically solve eq. (|19|) . The result is that V = is 
reached before all the local energy is gone only in the inner part of the V < region. The 
outer part cools completely while Yl is still lower than Y®, as it is apparent from Figs. 01 
There, the light thick lines represent the asymptotic lepton fraction and temperature 
profiles that would be reached if the antineutrino escape phase continued indefinitely 
without diffusion. In the r > 9 km region, the asymptotic temperature profile reaches 
zero. Correspondingly, the maximum value of Yl attainable before complete cooling is 
not Y® anymore. In practice, complete cooling does not occur. The profiles reached after 
one second, before diffusion becomes essential, are shown by the thin lines in the V < 
regions of figs. 01 The lines correspond to m 2 R = 10~( 10,9 ' 8 ' 7 ' 6,5 ) eV. The three largest 
values, that are not in the parameter space we are interested in, illustrate the asymptotic 
behavior of the profiles. Alternatively, they can be considered as the profiles at times 
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Figure 3: Lepton fraction and temperature profiles in the initial phase of the evolution. 
The thick dashed lines (dark solid in the V < region) represent the initial profiles. The 
dark thick lines show the profiles at the end of the neutrino escape phase. The thin lines 
show the effect of antineutrino escape 1 second later for m 2 R = l0 - ( 10 ' 9 ' 8 ' 7 ' 6 ' 5 ) eV. They 
can also be considered as the profiles at times t = lO"^ 5 ' 4 ' 3 ' 2,1 ' ^ sec (m 2 i?/(10 -5 eV)) 
provided that t < lsec. Finally, the light thick lines show the values that would be asymp- 
totically approached if diffusion could be indefinitely neglected. 



t = 10 - ( 5 ' 4 ' 3 ' 2 ' 1 ' °) sec (m 2 R/ (10~ 5 eV)) for a given value of m 2 R — as long as t < 1 sec, of 
course. The V = condition is reached and antineutrinos are locked in a portion of the 
initially V < region which is larger the larger is m 2 R. In a sufficiently large but finite 
time that would happen in all the region where the asymptotic Yl profile is constant and 
the T one is non vanishing. In the outer part of the inner core, on the other hand, neither 
Y® nor the asymptotic profile will ever be reached. However, the attraction of the fixed 
point on Yl has the effect of raising the Yl profile, with interesting implications that will 
be discussed in section EP1 The total energy loss in this phase is again a small fraction of 
the total available energy. 

After the first, fast phase of neutrino escape and the second slower phase of antineutrino 
escape, diffusion begins to play a significant role. This happens on the slowest time scale 
tdis- Diffusion "unlocks" the energy stored in the V = region. Some of it gets lost in 
the bulk, the rest diffuses out in the outer part and eventually is emitted mainly as active 
neutrinos. The diffusion regime however is different from the standard one that takes place 
in absence of transitions into the bulk. We have in this case a mixed regime in which the 
evolution is determined by an interplay of diffusion and escape. Let us consider first the 
inner part of the core, where V = 0. Due to diffusion, the densities change, thus spoiling 
the V = condition. The lepton fraction decreases and V becomes slightly negative. 
As V becomes negative, though, antineutrinos start to escape in the bulk thus giving a 
positive contribution to V . At some point, the positive contribution to V from conversion 
into bulk neutrinos will balance the negative contribution from diffusion. In this regime, 
some of the energy is lost in the bulk and some is diffused out. The amount of energy lost 
in the invisible channel turns out to be independent of m 2 R as long as the conversion of 
antineutrinos (slower than the neutrino conversion) is efficient enough to keep |V| small. 
The evolution in such a regime is well described by a single equation for the temperature. 
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In fact, denoting Yl = Y® + SYl, we have 



V = V(p,T)~^-5Y L . (20) 

When inserted in eqs. (|15jl. the previous expression allows to recover 5Yi(p,T) from 
eq. p5h). Eq. Q15b) then becomes 



n B T^ = a e (v ^ +V • Mo e + + a T )^-VT 2 



+ Ut„ - oT— — — - V • a e V/x^ - n s 



F 3/2 (ap u /T)J V V " / ^ 



(21) 



Eq. (|21|) always holds. In the small \V\ limit, one can approximate Yl ~ Y® everywhere, 
in which case the equation becomes selfconsistent and the evolution becomes independent 
of m 2 R. This of course will not hold forever since the efficiency of antineutrino conversion 
decreases as the star cools. Moreover, it certainly does not hold in the outer region where 
V is well below zero in the first place. 

In fig. |I] we show the evolution of the Yl and T profiles that follows from eqs. (|15[) 
for three values of m 2 R: lCT 10 eV, 1CT 9 eV and lCT 8 eV. The case of no new physics is 
also showed in each plot for comparison (dashed lines) . We do not aim at reproducing the 
results of the most sophisticated analysis in the literature for the latter case. However, a 
comparison between the solid and dashed lines in figure well illustrates the effect of new 
physics on the evolution. Due to neutrino conversion, Yl quickly drops to Y® in the inner 
part of the core. Then, on a longer time scale, diffusion further lowers Yl, thus starting 
antineutrino conversion. For large values of m 2 R, the latter is efficient enough to keep Yl 
close to Y® for a few seconds, until the region becomes too cool to sustain the necessary 
conversion. After the initial dramatic fall, the lepton fraction is higher than in the case of 
pure diffusion. The same happens in the outer part of the core, where the effect is more 
pronounced because neutrino loss never took place and \V\ is larger during antineutrino 
conversion. The lepton fraction even grows in the first second and then decreases slower 
than in the case of pure diffusion, at the expenses of a quicker cooling. The consequences 
of this effect will be discussed in the next subsection. For m?R = 10~ 8 eV, the temperature 
in the inner core essentially never increases. 



3.4 Implications for supernova physics 

Fig. 0] shows that deleptonization and cooling do take place on the same time scale as in 
absence of new physics despite the time scale of neutrino disappearance can be orders of 
magnitude faster than the diffusion one. What cannot be inferred from fig. |1] is the size 
of the portion of the available lepton number and energy that disappears in the bulk. Is 
something left to give rise to the SN 1987A neutrino signal and to revive the shock in the 
delayed explosion scenario? This issue is addressed in fig. El where the energy and lepton 
number lost by the inner core in the first 10 seconds are plotted against m 2 R and split 
in the component that goes into the observable neutrino flux and the component that is 
irremediably lost in the bulk. The latter grows with m 2 R but is never larger than 25% 
of the total. The total energy loss also grows with m 2 R. This is just a consequence of 
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Figure 4: Lepton fraction and temperature profiles at t = (0.01, 1, 10, 20, 40) sec for dif- 
ferent values of m 2 R (solid lines). The profiles for m 2 R = are also shown for comparison 
(dashed lines). 
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Figure 5: Energy and lepton number loss in the first 10 seconds. The total loss (solid 
lines) is shown together with its visible (upper dashed lines) and invisible (lower dashed 
lines) components as a function of m 2 R. 

the faster evolution of the temperature profile due to the additional loss (see fig. HJ). The 
net effect of the growth of both the total and the invisible energy loss is that the energy 
emitted in the visible channel in the first 10 seconds (as well as in the following 10 seconds) 
remains within 10% of the m 2 R = value. 

The lepton number lost in the bulk in the first 10 seconds is negative, a consequence of 
the dominance of the antineutrino escape phase. The reluctance of the lepton fraction to 
get away from the V = value shows in fig. El in the reduction of the total lepton number 
loss with m 2 R. It is also apparent how the large negative loss in the invisible channel 
translates into an enhancement up to a factor four of the emitted lepton number. Fig. 
shows that the enhancement takes place on the time scale of the initial antineutrino loss 
rate and is therefore effective already in the first second. 

The enhancement of the emitted lepton number has interesting consequences on the 
SN phenomenology. In order to discuss them, we first comment on the temperature of 
the neutrinosphere and the spectrum of the neutrino flux. A quantitative analysis of 
this and the following issues would require solving the evolution of the mantle, which is 
beyond the scope of this paper. We will therefore confine ourselves to some qualitative 
considerations. The disappearance probability is negligible at the density and temperature 
T v of the neutrinosphere, thus the latter can hardly be affected. Moreover, T u might even 
increase due to compressional heating until the mantle settles. The faster inner core cooling 
will therefore affect T v only later on. Finally, as the temperature of the neutrinosphere 
decreases, the neutrino mean free path grows and the neutrinosphere moves in the hotter 
interior, partially compensating the temperature change. We do not expect the spectrum 
of the neutrino flux to be significantly affected either. 

The composition of the neutrino flux depends in first approximation on T v and on the 
energy radiated per unit of lepton number. Neglecting the difference among the neutrino 
spectra |25| . the number of neutrinos emitted per unit of lepton number is N ~ kEj (LT U ), 
where E is the energy radiated with the lepton number L and k characterizes the shape 
of the energy spectrum at the neutrinosphere. In absence of new physics one has N ~ 10. 
Since the total electron lepton number of those neutrinos must be 1, there is a slight 
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number flux over the m 2 R = case. 
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Figure 7: Limits on the amplitude of v e , 
v e and Va T or &ut oscillations into an in- 
visible channel as functions of the squared 
mass difference Am 2 . 



prevalence of electron neutrinos over the other species. Consider now our situation in the 
extreme limit m 2 R = 10 -8 eV. The lepton number emitted in say the first second is 4-5 
times larger than in the standard case, while the energy emitted is about the same. Each 
electron neutrino emitted will be therefore accompanied in average by only a couple of 
neutrinos of other species. Such a peculiar neutrino flux composition could be indirectly 
observed in the existing neutrino detectors when the signal of a SN explosion in our galaxy 
will reach us. 

The prevalence of electron neutrinos in the neutrino flux for large m 2 R has two further 
consequences. First, it reduces the electron antineutrino component, which is constrained 
by the SN 1987A signal. Moreover, the electron neutrinos are much more effective than 
muon or tau neutrinos in depositing energy in the matter outside the neutrinosphere. As a 
consequence, a larger v e component might help a stalling shock in ejecting the SN envelope. 
In fact, the heating rate of matter irradiated by a hotter v e flux is proportional to T 2 , 
which we expect not to be significantly affected by a large m 2 R, and to the u e luminosity. 
As the left panel of fig. El shows, the total luminosity in all neutrino species (integrated 
over 10 seconds, but the same holds for the luminosity itself) is not significantly affected. 
Therefore, an increase of the v e component has the effect of increasing the v e luminosity 
and, in turn, the heating rate. The enhancement can be significant and is only limited by 
the reduction of the v e flux and the effect on r-processes one is willing to accept. 

4 "Conventional" invisible channels 

Before concluding, we discuss the limits on conventional (non self-limiting) invisible cooling 
channels. This simple application of the formalism set up in the previous section is relevant 
to our discussion because it provides the condition under which subleading contribution 
to the disappearance probability are under control. 

Let P U (L,E) be the probability that a neutrino of type v = u e , tv, u T , i> e , P„, 9 T with 
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energy E disappears in an invisible channel while traveling a distance L and let us as- 
sume that there is no reappearance 12 . In order to obtain a limit on P v , we compare the 
consequent lepton number loss rate with the diffusion one. The simplest case is that of 
a power-law disappearance probability P y (L,E) = Pq(L/ Lo) a (E/ Eq)P . For example, the 
mixing with a tower of KK states gives rise to a contribution to the escape probability 
which is essentially constant in the relevant range for E\ Pv e ,v e ~ (mR) 2 . 
The neutrino number density loss rate is 



(P)e 



(2vr)3 \ V (E) 



fv(E) , (22) 



where (P) = J °° dxP(x\ u (E), E)e~ x is the average over the distance traveled of the disap- 
pearance probability. In the case of a power law probability and within the approximations 
discussed in Section l.-S. 1| one has 

T » = T^ T5F4+ P- 2a ( [ T L ) ' (23) 

where a u = \ V E 2 /(6tt 2 ). By definition of idiff, the diffusion rate is T^is = nsdYL/dt = 
ubYl /tdiff- The limit on the probability following from T v < F^is reflects the abundance 
of the neutrino type considered, 

/ \ / \ \ 11 /200MeV\ 5 p Y L lOsec p . t . 

V Pv J Po 0.3 t diS 

Pv{\[T),T) < 1.5 • 10- V ( ±Wt^ fOT " = OT ^ ( 24b ) 

V 1 J PO 0.3 Ediflf 

P„(A(r),T)<3^0-X j ^- 2 »' 3 )(^) 5 ^|i^ = (24c) 

V T ) po 0.3 idiff 

where a a p = (5 + p - 2a)/(5T(a + 1)), b a/3 = 24/(T(5 + 0- 2a)T(a + 1)) are nu- 
merical coefficients normalized to aoo = &oo = 1- in particular, in the case of a con- 
stant electron neutrino disappearance probability P Ue = (mR) 2 we get the limit mR < 
0.5 • 10- 5 (200MeV//i l/ ) 5 / 2 (10sec/t d iff) that we used to determine the conservative upper 
bound m 2 R < 10 -8 eV. 

For completeness, we also give the limit on the amplitude of a simple oscillation prob- 
ability P U (L,E) = ylsin 2 (Am 2 L/(4£')) into an invisible channel. The limit on yl as a 
function of Am 2 for the case of u e , v e and v„ T or P^ jT oscillations is plotted in Fig. for 
p = po , Y = 0.3, t diS = 10 sec, p v = 200 MeV and T = 30MeV. 



5 Summary 

The supernova plays an important role in the phenomenology of a possible mixing between 
SM neutrinos and fermions propagating in the bulk of large extra dimensions. We studied 
both the constraints on such a mixing following from the necessity of avoiding unacceptable 

12 In the case of oscillations, reappearance can be neglected when P<1. The case of large mixing is 
not covered here. 
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energy loss and the implications for supernova physics. We considered the simple case of 
mixing between the electron neutrino and a single KK tower of 5D bulk neutrinos. The 
parameter space where interesting effects arise is quite broad, with the size of the extra 
dimension in the range 10 -2 eV < < IkeV and the mixing mass parameter m such 
that 1CP 12 eV < m 2 R < 10 -8 eV and mR < 1CP 5 . The interplay of neutrino diffusion and 
conversion into bulk neutrinos has an important role in the phenomenology we study and 
requires taking into account the relevant aspects of neutrino transport, which we did in 
the context of a simplified model of the protoneutron star core. 

Due to the large number of available states and especially to matter effects, the rate 
at which neutrino disappearance in the bulk cools the protoneutron star is initially dan- 
gerously high. If taken at face value, such a rate would set quite a stringent bound on 
m 2 R, m 2 R < I0 -12 eV in the case of neutrino conversion. However, the disappearance 
rate quickly reduces itself to acceptable values. This happens for different reasons. In 
the region where the MSW potential V is positive, neutrinos escape in the bulk and the 
lepton fraction quickly drops to the value Y® ~ 0.3 at which V = 0, thus stopping the po- 
tentially most dangerous MSW enhanced conversion before a significant amount of energy 
is lost. Then, on a much slower time scale, neutrino diffusion spoils the V = condition 
and unlocks the energy frozen in what was initially the V > region. The potential be- 
comes slightly negative and antineutrinos start escaping in the bulk to restore the V = 
condition. This feedback keeps Yl close to Y® until the temperature becomes too low 
to sustain the necessary antineutrino escape rate. In the region where the potential was 
initially negative, neutrino conversion never takes place. The initial loss rate is smaller 
because of the much lower antineutrino density. Still, it would be too high for large values 
of m 2 R. However, in the inner part of the V < region the antineutrino loss stops itself 
by drawing V up to zero before all the local energy is lost, analogously to the V > 
case. In the outer part of the V < region, there is not enough energy to reach V = 0. 
However, the conversion again controls itself. This happens this time because the energy 
loss, proportional to T 7//2 , becomes less and less effective as the region cools. 

The described mechanism allows to gain four orders of magnitude in the allowed range 
for m 2 R (10 -12 eV — ► 10 _8 eV). In the subsequent four orders of magnitudes, the MSW 
enhanced conversion would still be under control but what were previously neglected as 
subleading contributions to the oscillation probability become too large. Such contri- 
butions, as well as the limits on other "conventional" invisible cooling channels, were 
discussed in Sectional In the gained portion of the parameter space, the implications for 
supernova physics are particularly interesting. Although compatible with the SN 1987A 
signal, the evolution of the protoneutron star is significantly affected, especially for large 
values of m 2 R. Deleptonization and cooling take place on the same time scale as in absence 
of new physics. However, the reluctance of the lepton fraction to get away from Yl ~ 0.3 
slows the deleptonization, while the energy loss accelerates the cooling. Up to a fourth 
of the energy lost by the protoneutron star goes to the bulk but the total luminosity in 
all neutrino species turns out to be approximately the same as in the standard case, at 
least in the first 20 seconds. Antineutrinos escape in the bulk in the effort of keeping 
Yl close to Y®, which boosts the lepton number radiated from the neutrinosphere. A 
quantitative study of the consequences of such an enhancement would require solving the 
evolution of the mantle. However, we expect it to result in an enhancement of the v e 
component of the neutrino flux and therefore of the v e luminosity, an effect observable in 
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the existing neutrino detectors. Since we expect the temperature of the neutrinosphere 
not to be significantly affected, the energy deposition in a possibly stalling shock would 
also be enhanced at the expenses of a smaller v e luminosity. 

The exotic phenomenology discussed in this paper represents an example of how the 
standard ideas about protoneutron star evolution could be affected in an unexpected way 
by the presence of new physics. 
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